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1. Why submillimetre surveys? 

In this paper we briefly describe a 'key' survey at submm wavelengths which 
we are currently conducting to address some of the most important questions in 
cosmology - how, at what epoch and over what period of time did massive galaxies 
form at high-redshift? 

The primary motivation for undertaking surveys of high-z galaxies at submm 
wavelengths is the expectation, which has still to be proved, that the most massive 
galaxies {Mharyonic > 5 X lO^^M©) form the majority of their stars in a relatively 
short (< \Gyr), but extremely luminous phase at rest-frame FIR wavelengths and 
hence, at redshifts z > 3, at wavelengths in the atmospheric windows between 
350/j, — 1300/^ which are accessible from the ground. 

Throughout this paper we are referring to the rest-frame FIR emission when 
discussing the observed submm spectral energy distribution, and that all physical 
quantities are calculated assuming Hq = 50kms~^Mpc~^ and go = 0.5. 

1.1. EVOLUTIONARY STATUS OF HIGH-Z GALAXIES 

Submm and mm continuum observations offer a unique opportunity to determine 
some fundamental properties that provide a measure of the evolutionary status of 
high-z galaxies, but only after making the following assumptions: 



(a) that emission at submm wavelengths is due to optically-thin thermal re-radiation 
from dust grains at temperatures oiTrest — 30 — 70K, and is neither the optically- 
thin emission of a lower- frequency radio synchrotron component [e.g B20902+34 
Downes et al. 1996, Yun & Scoville - priv. comm.), or the optically-thick self- 
absorbed emission from a higher-frequency mid-IR synchrotron component {e.g. 
de Kool & Begelman 1989, Schlickeiser et al. 1991); 
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Figure 1 . The physical properties of high-z radio galaxies (open stars) and radio-quiet quasars 
(filled circles), which lie within the parallelogram, have been derived from submm or mm contin- 
uum detections (assuming T = 50K, M/^^/^^d = 500) and are compared to those of star burst 
galaxies (filled stars), ULIRGs (crossed circles) and elliptical galaxies (filled squares) in the local 
universe. The diagonal lines indicate constant hpju/yiu^. The vertical dashed-line shows the 
gas mass boundary, to the right of which is a region of the parameter space marked "PGs" where 
one can expect to find the projenitors of the most massive, > 5 X 10^^ Mq, elliptical galaxies. 
This figure is taken from Hughes, Dunlop & Rawlings (1997). 



(b) that the thermal emission from dust grains at an observed temperature Tabs is 
due to a high-z target, where Tabs = Trest/{^ + z), and is not confused with colder 
foreground galactic cirrus at 15 - 25 K; 

(c) that within the high-z galaxy, the grains are heated by young, massive stars in 
active starforming regions, and not by emission from the associated AGN. 

Given the above, it is possible to use the measured submm flux densities from 
high-z galaxies to estimate their rest-frame FIR luminosities, which are propor- 
tional to the starformation rate of massive stars at some early epoch. Additionally 
the remaining molecular gas mass available for future star-formation can be de- 
termined from a measure of the total dust mass, and compared to the expected 
baryonic mass in a present-day counterpart. For example observational evidence 
suggests that the most luminous host galaxies (L > 5L^) of low-z quasars and 
radio galaxies have masses > 5 x 10^^ Mq (Taylor et al. 1996), and thus we can 
describe the progenitors of similar galaxies forming around luminous high-z AGN 
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as primcBval if they can be shown to contain a high gas content, comparable to the 
stellar mass in low-z quasars and radio galaxies, and to be undergoing an intense 
burst of starformation. The uncertainties in the dust temperature, grain models, 
gas/dust ratio and their consequences in the determination of the FIR luminosities, 
SFRs, and gas masses in high-z galaxies have been described in detail elsewhere 
by Hughes (1996) and Hughes, Dunlop & Rawlings (1997). 

A summary of the physical properties of the few high-z radio galaxies and 
RQQs detected at submm and/or mm wavelengths is presented in fig. 1. These 
properties arc compared to those of the most himinous and active galaxies in the 
local universe. A conservative conclusion is that it is difficult, except possibly in 
the case of BR1202— 0725 (although it may be lensed) to describe any of the high- 
z galaxies detected at submm and mm wavelengths as gcinuincly prim,(Eval, since 
they lie outside the region of parameter space where, according to our definition 
of M_f/2 > 5 X 10^^ Mq and a SFR> 5OOM0yr"\ young galaxies should lie. Vary- 
ing the assumed dust tcnipcratiu'c through a reasonable range {2>QK 70 K), 
as illustrated by the representative locus in fig. 1, struggles to change this basic 
conclusion, but it does demonstrate the need to constrain the dust temperature 
which can only be achieved by making submm surveys at two or more wavelengths 
that straddle the rest-frame FIR peak at the highest redshifts (see §1.2). 

1.2. TECHNICAL FEASIBILITY OF FUTURE SUBMM SURVEYS 

It's all very well describing the importance of submm observations of high-z sources, 
but do the sensitivities of current and future instruments and telescopes make the 
detections of large numbers of high-z galaxies a realistic possibility? Such concerns 
are well justified given the extreme diffculty found in detecting only a handful of 
high-z AGN at submm wavelengths without the assistence of gravitational lensing 
(Dunlop et al. 1994, Ivison 1995, Isaak et al. 1994, Chini & Kriigel 1994, McMahon 
et al. 1994, Cimatti & Freudling 1995, Hughes, Dunlop & Rawlings 1997, Omont 
et al. 1996). 

In fig. 2 we illustrate the dependence of flux density on redshift for a typi- 
cal starburst galaxy SED in the major submm and mm atmospheric windows. 
The particular example shown is that of Arp 220 (log hpm ^ 12.3 Lq, SFR 
^ lOOM0yr"^), although the flux densities can be scaled for a starburst galaxy 
with any arbitrary luminosity or SFR at any redshift. The; conclusion that can be 
drawn from fig. 2 is that submm and mm instrumentation must achieve NEFD's 
of < 30mJyHz-i/2, < SOmJyHz-i/^, < 250mJyHz-i/2 at 1.3mm, 850/xm and 
450/im respectively if future studies are to detect, with reasonable integration 
times, high-z galaxies undergoing starformation at rates of > lOOMoyr"-'^. Such 
high and sustained SFRs are necessary if a massive galaxy, ^ lO^^Mo, is to be 
built in < IGyr - a timescale which, at high-z, corresponds to a redshift interval 
Az 4.5 — 2.5. These necessary instrumental sensitivities can be compared with 
those in table 1 for all the major existing (and proposed) submm telescopes. It is 
immediately clear that only the JCMT at 850/im and IRAM at 1.3mm currently 
provide sufficient sensitivity to carry out extensive cosmological studies, and that 
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further advances must await the future developments of the SMA, the South Pole 
10-m project, FIRST and the ambitious large mm arrays {e.g. MMA and LMSA). 

In addition to the gains in sensitivity, the availabiUty of new instrumentation 
on ground-based telescopes (e.g. JCMT, IRAM) and satellites (ISO) gives rise to 
a combination of improved resolution, imaging capability and greater wavelength 
coverage in the FIR, siibmm and mm. Therefore it should soon be possible for the 
first time to properly test the validity of the assumptions outlined in §1.1, partic- 
ularly quantifying the level of cirrus confusion on scales of ~ 10", constraining the 
dust temperature, and discriminating between the competing thermal and non- 
thermal emission mechanisms in the FIR by measuring directly, and with high 
photometric precision, the rest-frame submm spectral indices of high-z galaxies 
(Chini et al. 1989, Hughes et al. 1993). Finally, since all massive galaxies at liigh-z 
have been originally identified by detecting a luminous AGN then the effect of the 
AGN continuum emission, through primary or secondary processes, on the FIR 
luminosity and on the overall evolution of the host galaxy must be quantified. 
This can be achieved by either detecting genuinely non-AGN galaxies, previously 
unidentified, in new submm blank field surveys, or by making pointed submm ob- 
servations of samples covering a broad range of AGN luminosity and redshift. We 
describe the latter of these alternatives below. 
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TABLE 1. Continuum sensitivities of receivers on submm/mm tele- 
scopes. The SMA(6) NEFDs assume 6 antennae in the array; SMA(8) 
includes two additional Taiwanese antennae, and in parentheses the use 
of dual polarization heterodyne receivers at 345GHz. Numbers in italics 
are predicted sensitivities. 







NEED (mjy Hz" 


-1/2) 


telescope 


area (m^) 


1350Mm 


850/im 


450/im 


JCMT 


177 


60 


80 


700 


CSO 


79 






2000 


IRAM 


707 


60 






SEST 


177 


200 






SOFIA 


5 




200 


200 


SMA(6) 


170 


102 


238 


2040 


SMA(8)+JCMT+CSO 


481 


36 


84 (59) 


719 


South Pole 10-m 


79 




50 


60 


MMA 


2010 


8 


25 


150 


FIRST 


7 




54 


54 



2. Submm observations of radio samples covering the P-z plane 

We have chosen to address some of the questions of evolution and formation of 
massive galaxies by observing complete sub-samples of radio galaxies that span a 
factor of ^ 1000 in radio luminosity (logP408MHz ^ 25.0 — 28.0 WHz~^sr~-'^) and 
range in redshift between 1 < z < 5. The samples drawn from 3C, 6C, 7C and 
LBDS surveys are shown in fig. 3, and taken together eliminate the tight correlation 
between radio luminosity and redshift observed in any single survey by covering 
the entire P-z plane. It is well known that powerful low-z radio sources reside in 
elliptical galaxies and, assuming the same to be true at high-z, we can therefore 
study the gas mass fraction and level of star formation in massive host galaxies 
as a function of AGN luminosity at a particular redshift, and follow the evolution 
of these quantities over a period of ~ 4Gyr {i.e. from ^; 5 — > 1). We will report 
shortly on the first JCMT observations of these samples at 850/um and 450/im 
using the new bolometer array SCUBA. 
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Figure 3. The radio-luminosity:redshift (P — z) plane of liigh-z galaxies selected from various 
radio surveys (3C - open circles; 6C - filled squares: 7C - crosses; LBDS - filled circles) to be 
observed at submm wavelengths. Using these samples we can quantify the contribution of an 
AGN to the rest-frame FIR luminosity, and trace the evolution of gas mass and star formation 
rate as function of redshift and radio luminosity. 
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